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CN \ ABSTRACT 

d ■ 

The effective content of cosmic rays (CR) in galaxy clusters remains elusive. The evidence of relativistic electrons (RE) in the 
, subset of clusters endowed with a radio halo remains hardly quantitative in the absence of robust estimates of the magnetic 
f****) ■ field B(r), derived from Faraday Rotation (FR) measurements. The content in relativistic protons (RP) requires a different 
£SJ ' approach, the only direct one residing in the detection of their collisional production of gamma rays (GR). 

, Based on the evidence of merging phenomena in clusters, theory predicts a large content of RP, whose energy density could 
t-H ■ be a large fraction of the thermal energy. This paper aims to estimate a maximum production of both secondary relativistic 
^ ' electrons, SRE. and GR from the RP that have supposedly accumulated throughout the entire history of a cluster. 

SRE and GR production is maximized when the RP and the thermal gas share the same radial profile. The production rate 
is normalized by adopting a reference value of 0.3 for £, the ratio of RP to thermal pressure. The SRE content which obtains, 
when constrained to reproduce the observed radio brightness profile, yields univocally B(r), if the presence of primary RE were 
negligible. 

' This procedure is applied to four radio-halo clusters (Coma, A2163, A2255, A2319). In these objects, the central value -Bo 
f"*^ , required is consistent with typical, albeit rather uncertain, values derived from FR, although for A2163 and A2319 no reliable 
FR estimates are available to strengthen this result. On the other hand, B(r) typically increases beyond the thermal core, a 
hardly acceptable condition. This problem is alleviated by assuming a mix of SRE and of "primary" RE (PRE), with the latter 
becoming the dominant component beyond the thermal core. These results suggest that in clusters without a radio halo detected 
Q ' so far a diffuse radio-emission should also be observable due to SRE alone, and therefore more centrally condensed, provided 
\ that £ is of the order of 0.3. To encourage deeper radio observations of such clusters, some examples were selected that seem 
. rather promising. Efforts in this direction, if accompanied by FR measurements, could provide highly significant constraints on 
the CR content in clusters, even before the future GLAST mission will have accomplished the hard task of detecting the GR. A 
complementary result concerns the excess far UV in the Coma cluster, that some authors have attributed to IC emission from 
SRE. It is shown that this hypothesis can be excluded, because it requires a RP energy content in excess of the thermal one. 
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1. Introduction 

A significant fraction of galaxy clusters are endowed with 
a diffuse, low-brightness radio-emission, which is clearly 
distinct from the extended radio components, if present, 
associated with member galaxies. The historical prototype 
(Willson 1970) is the radio source Coma C in the Coma 
cluster. These so-called radio halos represent direct evi- 
dence of relativistic electrons and magnetic fields embed- 
ded in the thermal plasma responsible for the X-ray emis- 
sion. It is very likely that relativistic protons are present 
as well. In the following these components will be labelled 
RE (electrons), RP (protons), CR (for cosmic rays, pro- 
tons and electrons), and B (the magnetic field and its 
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strength). Since the brightness of the radio emission de- 
pends (non-lincarly) on B and the RE density, it cannot 
be excluded that the presence of a non-thermal compo- 
nent is common to all clusters, including those without, 
as yet, an observed radio-halo. 

Generally speaking, we are still far from solid quan- 
titative estimates of the non-thermal content. The radio 
emission alone can only provide a minimum for the energy 
content in CR and B combined, according to the recipe 
first proposed by Burbidge (1956). This minimum is well 
below the thermal energy content (Govoni & Feretti 2004). 

A direct way to estimate the RE and B contents sep- 
arately was first applied by Perola & Reinhardt (1972) to 
Coma C. Although at the time they could (incorrectly, as 
it became clear later) and did attribute the whole of the 
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X-ray emission to this process, the RE inevitably lose en- 
ergy on the photons of the cosmic microwave background 
through the so-called inverse Compton (IC) process, thus 
it remains obviously true that, within a rather wide range 
of B values, the energy lost should be observed in the X- 
ray band. All the attempts, made so far, to disentangle this 
contribution from the thermal X-ray emission, which have 
concentrated attention on the hard (greater than 10 — 30 
keV) regime, have led to results that are rather contro- 
versial particularly in their interpretation (Fusco-Femiano 
et al. 2004; Rossetti & Molendi 2004). This follows from 
the evidence of an excess, possibly IC emission, obtained 
with non-imaging instruments, which therefore can only 
be used to estimate the volume-integrated amount of RE 
within the instrument's field of view and, rigorously speak- 
ing, only in a very narrow energy interval, 1.6 — 3.2 GeV. 
In a previous paper (where most references can be found, 
Colafrancesco, Marchegiani & Perola 2005, Pap. I), we 
discussed this issue at length and emphasized the need 
to measure the surface brightness distribution of the IC 
emission, in order, first, to verify whether the RE respon- 
sible for this radiation are effectively confined within the 
observed radio halo and, second, to solve the degeneracy 
intrinsic to the synchrotron emission and therefore obtain 
fundamental information on the radial dependence of the 
RE density and of B. 

The origin of the RE was outside the goal of Paper 
I. They can have a "primary" origin, which is accelerated 
along with the RP. However, there must also be RE of 
"secondary" origin, produced through the inevitable in- 
elastic interactions of the RP with the ambient plasma. 
This process, notably, gives rise to the emission of gamma 
rays (GR) with energy > 50 McV. If we knew the number 
and distribution of the RE (from their IC emission), the 
question (Dennison 1980) whether the "secondary" ones 
(SRE) are either a negligible or an important fraction of 
the whole RE population would require assumptions on 
the amount of RP, with basically one single important 
constraint, namely that their pressure should not exceed 
that of the thermal gas. Vice versa, if by measuring the 
spectrum, flux, and the associated surface brightness of 
the gamma rays (after verifying that they cannot be due 
to other mechanisms, such as the hypothetical annihila- 
tion process of the dark matter, see for instance Pieri & 
Branchini 2004, Colafrancesco et al. 2006) one could derive 
the amount of RP, a very stringent constraint on the pro- 
duction rate of the SRE could be placed. If this rate were 
too low for a significant contribution to the radio-emission, 
except for unacceptably high values of B, the answer to 
the question would be immediate. It is clear, however, that 
only the combination of reliable measurements of the GR 
and the IC X-rays, with the latter leading to an absolute 
estimate of the RE amount and distribution, will close 
the circle. The amount and distribution of the RE could 
also be obtained via an independent estimate of B (radial 
distribution included), which can be derived from Faraday 
rotation (FR) measurements, once this goal is reached in 
an uncontroversial way (see Govoni & Feretti 2004 for a 



discussion of the state of this art). In either (IC or FR) 
case, the GR represent an irreplaceable step for the issue 
of the CR content. 

In this context, the contribution that the present paper 
aims to provide, given the still immature observational 
status, consists in evaluating the "maximal" requirements 
on the RP and on B, including their radial distribution, 
assuming that the radio-emission is "totally" due to the 
SRE. These requirements, constrained to avoid conflicts 
with either physical arguments (the non-thermal pressure 
should not exceed the thermal pressure) or observational 
evidence (limits on B from FR), will be used to predict 
the uppermost value of the GR emission intensity, which 
one would expect to inevitably follow. 

The authors are aware of the fact that arguments have 
been put forward that are against a dominant contribution 
to the RE of the SRE component (e. g. Brunctti 2003). 
Nonetheless, the approach adopted is useful on two dif- 
ferent grounds. First, the RE lifetime in the energy range 
of interest is short (about one tenth) compared to the 
age of a cluster; thus the currently relevant primary RE 
(PRE) must have been accelerated in this last fraction 
of the cluster age, while the RP giving rise to the SRE 
in the same fraction of time have likely accumulated over 
the entire cluster lifetime, meaning that their contribu- 
tion to the RE content might be far from negligible with 
respect to that due to the acceleration processes. Second, 
the choice of maximizing this contribution leads in a sense 
to the most optimistic predictions on the intensity of the 
GR emission, which could be immediately tested by its 
measurement. 

Another item is connected with the excess emission re- 
ported in the far ultraviolet-soft X-rays in the Coma clus- 
ter in particular, that might be attributed to IC emission. 
If this were the case, the relevant RE population would 
fall into a totally different range of energies with respect to 
the radio-emitting ones. In this hypotheses, some authors 
(Bowyer et al. 2004) have attributed a "secondary" origin 
to these RE. Its relevance to the ensuing requirements on 
the parent RP will also be discussed. 

The paper is organized as follows. Section 2 contains 
a brief overview of the theoretical results on the acceler- 
ation of RP in clusters. Section 3 is devoted to the as- 
sumptions and prescriptions on which the subsequent cal- 
culations are based. Section 4 illustrates the properties of 
the four clusters used as examples and the way their radio 
brightness distribution is described analytically. In Section 
5 the results are reported for the GR emission and for the 
strength and radial distribution of B, when the radio halo 
is "totally" attributed to the SRE. Section 6 is devoted 
to the issue of the UV excess in the Coma cluster as due 
to IC from SRE. Section 7 is devoted to a discussion, and 
Section 8 to the conclusions. 

2. Acceleration and accumulation of RP 

Since Jaffe (1977) discussed the possible diffuse origin 
(that is, in situ acceleration versus injection by radio 
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sources) of the RE in Coma C, science has progressed a 
long way. In the theoretical scenarios of structure forma- 
tion in the Universe through hierarchical clustering, with 
the relatively recent support from Chandra and XMM- 
Newton observations, the merging processes are widely 
suspected to be the events that can lead to a diffuse ac- 
celeration of CR in clusters, along with an increase in the 
energy content of the thermal gas. 

The papers quoted below do refer to the results 
(Blandford & Eichler 1987) of the shock diffusive accel- 
eration theory that predict spectral indexes s linked to 
the Mach number M as 



M 1 + 1 
! M 2 -1' 



(1) 



From hydrodynamical simulations (e. g. Ryu et al. 
2003) it appears that in the central regions, where the 
largest amount of particles are processed, the typical val- 
ues of the Mach number are 2 < M < 4, and, according 
to Eq. ([1]), those of spectral index are 2.3 < s < 3.3. The 
amount of energy going into the CR channel can be as 
large as 30 — 50% of what goes into the thermal chan- 
nel (Miniati et al. 2001a; Ryu et al. 2003). Much higher 
values of M, which are even more efficient in the acceler- 
ation process, are attained in the cluster outskirts, where 
(Pfrommer et al. 2006) a higher ratio between the two en- 
ergy contents could eventually prevail, except that dynam- 
ical stability problems might then ensue. In this paper, as 
explained in the next section, care is taken that nowhere 
within the radio halos, as defined by the outermost con- 
tour of the available maps, does the non-thermal pres- 
sure (inclusive of CR and B) exceed the thermal pressure. 
Furthermore, a single value of s is adopted for simplicity, 
namely s — 2.7. (Incidentally, this is the value required 
for reproducing with a SRE population the radio slope 
of Coma C, the only one well determined among the four 
cases which will be considered in the following). A value 
of s higher than 2 implies that the choice of the minimum 
energy is critical to the CR pressure. Following Gieseler, 
Jones & Kang (2000), the minimum adopted Lorentz fac- 
tor is given by 



7« 



1 



3.4 x 10~ 5 fcT, 



keV 



(2) 



where kTk e v is the thermal gas temperature, and C2 is 
approximately equal to 4. Thus jmin is practically equal 
to 1 in all cases. 

For the RP momentum (p) spectrum, the form 
adopted is N R p(p) oc p~ s , which can be written as a func- 
tion of 7: 



N RP {~/,r) = K RP {r)- 



7 



( 7 2 _ l)(s+l)/2' 

The associated pressure is given by 



(3) 



1 fv™» 

P RP (r) = -m p c J N RP ('Y,r)P'Yv('y)d'Y. (4) 



The radial dependence of the RP content, with a constant 
s, as anticipated, will be described through that of their 
Prp- 



3. Assumptions and prescriptions 

The production rates per unit volume of the SRE, 
Qe(le, r) and of the GR, Q g (E g , r), which are strictly re- 
lated to each other, are proportional to the product of the 
thermal and non-thermal proton density. Since the first 
of the two quantities is known, the radial behavior of the 
second, which maximizes the production rates, must be 
found. 

According to what has been described in Sect. 2, and 
assuming that among the CR the RP are the energeti- 
cally dominant component, it is convenient to introduce 
the ratio between P R p and Pth, the pressures of the non 
thermal and the thermal components, 



Prp(t) 



Ptn{r) ' 

which is bound to remain everywhere less than unity: 
£0) < 1. 



(5) 



(6) 



So long as kT can be regarded as constant (a reason- 
ably good description of clusters within the extent of their 
radio halos), and under the assumption that the mean RP 
energy is also constant (a direct consequence of the as- 
sumption of a constant value of the exponent s anticipated 
in Sect. 2), £(r) oc n R p/n t h, where n R p (= J Nppdj) and 
n t h are the number densities of the relativistic and the 
thermal protons, respectively. Without loss of generality, 
the radial dependence of n R p can be chosen with the same 
functional form as is generally used for nth, 



n t h{r) = nthfi 



r c .th 



-9th 



(7) 



with the two parameters, core radius and exponent, left 
free: 



n RP (r) = n RP<0 



r c , R p 



2~| ~1RP 



(8) 



Not surprisingly, the production rates reach their max- 
imum values when r CiR p — r c ,th and q R p — qth, therefore 
for £ constant with radius, which will be labelled £o- These 
values are given by 



Qe(le,r) = n tht0 K R Pfi 



(7 2 _ 1 j ) (s+ i)/ 2 ^ ± (7 P )^e(7e,7 P )rf7 P (9) 



r c .th 



1 —2qth 



c/3(7 P ) 
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Qg(E g ,r) = n t h,oK R Pfi 



( 7 a _ 1 )( s +i)/2 ^°(>) J? g( £; g'>) rf > (10) 



1 + 



r c ,th 



- —2qth 



where a v ± and a^o are the cross sections for the produc- 
tion of charged and neutral pions, respectively, in proton- 
proton collisions, and the functions F e and F g are given 

by 



^77 r ±^e(7e,77 r ±)^7r±(7 7 r±,7p) (H) 



where F e (pf e , 7^±) is the electron energy distribution func- 
tion from pions, and F v ± (7,r± , J P ) is the charged pion en- 
ergy distribution function; 



F g (E g n P ) 



d7 Jr oF s (£; s ,7 OT o)F 7r o(7 7r o,7p) (12) 



where F g (E g , 7^0 ) is the photon energy distribution func- 
tion from pions and F v o (7^0 ,j p ) is the neutral pion energy 
distribution function. The distribution functions in (jlip 
and (fl"2)) used here are taken from Moskalenko & Strong 
(1998). 

Given a value of £0, the GR spectral brightness and 
luminosity can be obtained immediately from integration 
of pop out to a maximum radius R. For the SRE, in order 
to obtain their number density as a function of energy and 
radius, n e (7 e , r), it is necessary to solve the equation that 
locally governs (spatial diffusion can be regarded as irrel- 
evant in this context, see e. g. Sarazin 1999) the balance 
between injection and energy losses, 



d 

— {b(-f e ,r)n e (j e ,r)] = -C 

<?7e 

whose solution is given by 



>Ale,r) 



n e ("f e 



QeK,T)di e . 



(13) 



(14) 



In (| 13[) and (fT4)) . b(~f e ) = —d, r y e /dt is proportional to 
the sum of the square of B mw b, the magnetic field corre- 
sponding to the energy density of the cosmic microwave 
background (IC losses) and the square of B(r). The latter 
quantity is not fixed a priori, but will be derived through 
an iterative process, when the SRE is requested to re- 
produce the radio brightness distribution for any given 
cluster. Then a radial behavior of the magnetic field will 
univocally emerge, and the normalization B(r = 0) = B 
will follow from the adopted value of £q- As a reference 
value for the maximum in the GR emission to be expected, 
£0 = 0.3 will be adopted as representative of the most ex- 
treme conditions described in Sect. 2 with regard to the 
RP acceleration. 



4. Four clusters and the brightness profile of their 
radio halos 

The exercise is applied to four clusters as examples. Their 
names and properties are summarized in Tables [1] and [2l 
In Table [I] in a sequence, after redshift and temperature, 
the three parameters n t h,o, fc,tht an d 1th are given, which 
quantitatively define the density distribution n t h, accord- 
ing to ((7|) . As in Paper I, the radial (azimuthally averaged) 
brightness distribution at frequency v was fitted with a 
functional form similar to the one usually adopted for the 
X-ray brightness, namely (the subscripts R here mean "ra- 
dio") 



S v {9) — S v .q 



'c.R 



-'In 



(15) 



where 9 is the angular distance from the center. Under the 
assumption of spherical symmetry, the profile of the radio 
emissivity is similar in form to n t h. ([7]), with qp — q' R + l/2 
and r c r — 9 C pD, D being the cluster angular distance. 
These are the two quantities measured at v = 1.4 GHz 
in Table [21 along with the flux at the same frequency, 
the radius Rh, and the angular size Oh of the radio halo. 
The spectral index ctp in the second column is fairly well- 
determined only for the Coma cluster; for the three other 
objects, the value is given with a double dot, to remind 
the reader that it has not been measured well yet (the flux 
density data available will appear in Figs. 1-4). Note that 
the value of s — 2.7, chosen in Sect. 2 as a reference energy 
spectral slope of the RP, at equilibrium yields a slope of 
the SRE that almost exactly matches the radio spectral 
index of Coma (FigJT]) . 

5. Numerical results 

The results are summarized in Table [3] and illustrated in 
Figs. [U [21 [3l [U and El The following quantities are given 
in Table [3J the number density of the RP at the center, 
nRPfi, the value of Bq for which the radio flux is best 
reproduced by the SRE alone, the gamma ray flux from ir° 
decay, / s (100 MeV), obtained integrating out to the radius 
of the radio-halo Rh-, the gamma ray luminosity L g in the 
band 0.1 to 10 GeV, and next to it, for comparison, the 
radio luminosity approximated as 1.4(GHz) xL(1.4 GHz). 
The quantity Bq in the fourth column is the reduced mean 
value of B if "small-scale" scalar fluctuations (see Paper I) 
are admitted to exist on the order of < {SB) 2 > / B 2 = 1 
throughout a cluster. 

In the upper panels of Figs.[T]to[H the radial behaviors 
of P t h, Prp, and of the magnetic field pressure are shown 
out to Rh ■ Note that the magnetic pressure remains lower 
than the other two everywhere. In the lower panels of the 
same figures, along with the observational data, the ra- 
dio spectral flux density is drawn for three values of Bq, 
the one given in Table [3] and two more, to illustrate the 
sensitivity of the normalization to this quantity. 

In Fig. [5] the spectral density distributions of the 
gamma rays from ir° decay of all four clusters are given. 
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kT 


n t h,o 






Qi h 






(keV) 


(hYo 2 cm- 3 ) 


(/ifo 1 kpc) 




Coma 


0.023 [1] 


8.2 [2] 


3.4 x 10~ 3 


[2] 


300 [2] 


1.125 [2] 


A2319 


0.0557 [1] 


9.12 [3] 


7.83 x 10" 3 


[4] 


94 [5] 


0.765 [5] 


A2255 


0.0806 [1] 


7.3 [6] 


1.72 x 10~ 3 


['] 


410 [8] 


1.11 [8] 


A2163 


0.203 [1] 


14.6 [9] 


7.87 x 10~ 3 


[9] 


220 [9] 


0.93 [9] 



References: [1] Struble & Rood (1999); [2] Briel, Henry & Bohringer (1992); [3] Arnaud & Evrard (1999); [4] Trevese, Cirimele 
& De Simone (2000); [5] Feretti, Giovannini & Bohringer (1997b); [6] David et al. (1993); [7] Jones & Forman (1984); [8] Feretti 
et al. (1997a); [9] Elbaz, Arnaud & Bohringer (1995). 

Table 2. The four clusters, parameters of the radio halo 





a R 


F(1A GHz) 


r c ,n 




Rh 


e h 






(Jy) 


{h^l kpc) 




(h^l kpc) 


(arcmin) 


Coma 


1.35 [1] 


0.64 [1] 


670 


4.4 


900 [1] 


64 


A2319 


1.8: [2] 


0.153 [2] 


37 


0.98 


470 [2] 


16 


A2255 


1.7: [3] 


0.043 [3] 


470 


3.0 


670 [3] 


16 


A2163 


1.18: [4] 


0.155 [5] 


220 


1.4 


1,050 [4] 


12 



References: [1] Deiss et al. (1997); [2] Feretti et al. (1997b); [3] Feretti et al. (1997a); [4] Feretti et al. (2004); [5] Feretti et al. 
(2001). 

Note that the expected flux scales linearly with £o, and 
the values in the figure correspond to £o = 0-3- For com- 
parison the sensitivity curves are shown for the experi- 
ment EGRET onboard the satellite Compton-GRO and 
for the future AGILE and GLAST missions, for an ex- 
posure time as reported in the caption. Although the 
curves apply to point-like sources, note that the angu- 
lar resolution (on-axis) of GLAST is about 60 arcmin at 
1 GeV (from http://www-glast.slac.stanford.edu), hence 
it should also apply to the extended sources considered 
here (see Table 2 for the angular size of the radio ha- 
los). It is evident that out of the four clusters considered, 
only the Coma cluster emission could be significantly de- 
tected by GLAST, even if the value of £o were lower than 
0.3, down to about 0.1. With £ = 0.3, Coma could be 
marginally detected by AGILE, while A2319 and A2163 
could be marginally detected by GLAST. 

The discussion of these results is postponed to Sect. 
7, after considering a case of emission, in the ultraviolet, 
which has been proposed in the literature as due to SRE. 

6. SRE and the EUV excess in the Coma cluster 

Within 25 arcmin of the center of the Coma cluster, the 
EUVE satellite detected an emission in the far UV band 
130-180 eV, as reported in Lieu et al. (1996) and Bowyer 
et al. (1999, 2004). It exceeds the extrapolation of the 
thermal emission measured in the X-rays. The origin of 
this excess is controversial (Lieu et al. 1996; Enfilin & 
Biermann 1998; Brunetti et al. 2001; Bowyer et al. 2004), 
here only the hypothesis put forward by Bowyer et al. 
(2004) (see also Kuo et al. 2005 for a generalization to 
the clusters as a class) will be considered, namely that, in 
the first place, the emission mechanism is IC on the mi- 



crowave background and, in the second, that the electrons 
have a secondary origin. The reason for inserting this is- 
sue into the present paper follows from the fact that the 
UV measurement, albeit essentially monochromatic and 
without spectral information, is spatially resolved; hence, 
under the IC and spherical symmetry assumptions, both 
the number and the radial distribution of the electrons, 
whose energy must be j e ~ 300, can be derived in ab- 
solute terms. The amount of parent RP can therefore be 
derived, and their pressure compared with the thermal 
pressure to verify if the secondary origin of the electrons 
remains tenable. 

The profile of the UV surface brightness Sjjv(9) can 
be described with the functional form of Eq. (fT"5|) . and it 
turns out (as pointed out by Bowyer et al. 2004) to be 
very similar to that of the thermal X-ray emission. A best 
fit with 9 Ct uv = 0c,x = 10.5 arcmin yields an exponent 
q' uv = 1.92 ±0.20 \x 2 red = 0.64 with 11 dof), consistent 
within la with q' x = 1.75 (Briel, Henry & Bohringer 1992). 
The emissivity is given (see Sect. 4) by a functional form 
like 0, with r CtUV = r c ,th and q uv = q' uv + 1/2 = 2.42 ± 
0.20: in the following qjjv = 2.42 will be used. Since the 
emissivity is proportional to the number density of the 
electrons n e with 7 e about 300, this quantity will have 
the same radial dependence: 



n e (r) = n e (0) 



r c ,th 



(16) 



The RE with j e about 300 would emit synchrotron ra- 
diation at v ~ i)A(B/ixG) MHz, an inaccessible window. 
As pointed out already by Bowyer & Berghofer (1998), 
the radial behavior (1161) is likely to be very dissimilar from 
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Table 3. Numerical results 



n RP: o B B / 9 (100 MeV) L g (0.1-10 GeV) vL R (1.4 GHz) 
(cm" 3 ) jfiG) (fiG) (cm" 2 s" 1 GeV 1 ) (erg s" 1 ) (erg s' 1 ) 



Coma 


1.1 x 10" 


-b 


1.2 


0.8 


8.4 x 10' s 


1.3 x 10" 


9.9 x 10 39 


A2319 


2.4 x 10" 


-5 


4.6 


3.2 


8.6 x 10 -9 


8.6 x 10 42 


1.5 x 10 40 


A2255 


6.0 x 10" 


(i 


4.0 


2.8 


2.2 x 10" 9 


5.4 x 10 42 


9.7 x 10 39 


A2163 


2.0 x 10" 


-5 


2.4 


1.7 


6.0 x 10 -9 


1.0 x 10 44 


2.6 x 10 41 



_ 



10° 



10 



1 10 



10 
10" 



-3 " 

-a : 

-5 : 
0.01 



0.10 
r (Mpc) 



E 
o 

> 





v (MHz) 

Fig. 1. Coma cluster, £o = 0.3, s = 2.7. Upper panel: Pth 
(full line), Prp (dashed line), and the magnetic field pressure 
(dotted line), with Bo = 1-2 /iG, as a function of distance from 
the center. The vertical dashed line marks the position of r c ^h- 
Lower panel: the SRE radio emission spectra labelled with 
three values Bo (in ^G). The data points are from Thierbach, 
Klein & Wielebinski (2003). 



that of the RE responsible for the observed radio emis- 
sion; otherwise, it would imply a magnetic field increasing 
with r out to at least 2/3 Rh for the radio brightness dis- 
tribution of Coma C. 

Given the SRE density distribution (fT6|) . one can re- 
sort to the population of the parent RP by inverting 
the procedure described in Sect. 4. The first step con- 
sists in obtaining Q e ("f e , r) from Eq. (|13p . Concerning the 




v (MHz) 

Fig. 2. A2319. Same as Fig. 1, except for: upper panel, Bo = 
4.6 fiG. Lower panel, data without error bars from Feretti et 
al. (1997b), with 20% error bars added arbitrarily. 



energy-loss rate, 6(7e, r), it now consists of two terms, the 
Coulomb loss rate 6 C , which depends on r since b c oc nth, 
and the IC loss rate bic, which is constant in space and 
corresponds to the loss rate in an equivalent magnetic 
field B mw ij = 3(1 + z) 2 /iG. The synchrotron losses de- 
pend on the unknown value and radial behavior of the 
magnetic field: ignoring these losses does not affect the 
conclusion, as will be clear later on. Using the expression 
for b c given in Sarazin (1999), it turns out that (for n t h 
constant with time) at r much less than r c and for elec- 
trons with 7 e ~ 300, b c was equal to bic at z* = 0.384, 
and their lifetime was 2.1 Gyr at that epoch. The time 



P. Marchegiani et al.: Cosmic rays in galaxy clusters 



7 




1 u [ 

0.01 0.10 

r (Mpc) 




Fig. 3. A2255. Same as Fig. 1, except for: upper panel, Bo = 4 
/J.G. Lower panel, data without error bars from Feretti et al. 
(1997a), with 20% error bars added arbitrarily. 
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Fig. 4. A2163. Same as Fig. 1, except for: upper panel, Bo = 
2.4 ^G. Lower panel, data from Feretti et al. (2001; 2004), 
whose error bars are within the height of the two triangles. 



elapsed since then up to the epoch corresponding to the 
redshift of the Coma Cluster - for Hq = 70 km s -1 Mpc -1 , 
O ro = 0.3, and Q A = 0.7 - is 3.4 Gyr, therefore the UV 
emission observed now from the innermost region of the 
cluster is due to electrons affected more by Coulomb than 
by IC losses. At larger distances from the center the latter 
dominate instead. 

Equation (fTB")) contains the spectral energy distribution 
of the electrons, which is unknown. The calculations were 
carried out in full, that is by using Eq. (|13[) , in combination 
with Eq.@, adopting two different values of s for the RP, 
namely 2.7 and 2.3. 

The radial distribution of the parent RP, which repro- 
duces the UV brightness profile fairly well (see Fig. [5]), is 
given by 



n RP {r) oc 



1+1- 



2~| ~QRP 



(17) 



with r c = r Ct th and qpp — 2 (irrespective of which one of 
the two values of s is used). The exponent qpp is rather 
well constrained: as shown in Fig. [5J the UV brightness 
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Fig. 5. Spectral energy distribution of the gamma ray flux 
from 7T° decay expected in the case illustrated for each clus- 
ter in Figs. 1 to 4. The sensitivity of EGRET, AGILE, 
and GLAST is drawn for comparison purposes and corre- 
sponds to a 5a detection in a full-year sky survey (from 
|http://people.roma2.infn.it/~ agile / ) . 
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Fig. 6. The UV brightness profile from SRE for two values of 
qnp (see text). The data are from Bowyer et al. (2004). 



Table 4. The value of £(r = 0), for two choices of the expo- 
nents s and qnp, needed to explain the UV excess in the Coma 
Cluster as due to SRE 



s — 


2.7 £(r = 0) 


s — 


2.3 £(r = 0) 


qnp = 


1.125 9 


qpp = 


1.125 7 


qpp 


= 2 21 


qpp 


= 2 15 



distribution could not be reproduced if the RP shared 
the same radial distribution with the thermal gas, namely 
if Qrp — Qth — 1-125. This would, notably, turn out to be 
the case if the Coulomb losses had been ignored and only 
the spatially "constant" IC losses taken into account. The 
result obtained implies that the quantity £(r), defined by 
([5]) is not constant and that it is maximum at the center. 

The normalization of Eq. (fl7|) and the corresponding 
Prp depends on the exponent s. The values of £(r = 0) 
for the two choices of s, obtained by integrating the RP 
spectrum down to the 7 m j n given by Eq. ^ , are given in 
Table [Hand are much higher than unity for both values of 
qnp in Fig. O The inevitable conclusion, which would be 
even stronger if the synchrotron losses were non-negligible, 
is that the hypothesis of a secondary origin of the electrons 
is not tenable in this case. 

Surprisingly, Bowyer et al. (2004) reach the opposite 
conclusion. This follows, partly, from their decision to ig- 
nore the Coulomb losses, in addition to the synchrotron 
losses. But this implies at most a difference in the value of 
£ by a factor five. Since other approximations only have 
marginal effects, one is led to conclude that their result 
is affected by a material error. In this respect, it is worth 
mentioning that Miniati et al. (2001b) conclude that the 
UV excess due to SRE should be lower than the measured 
value, on the basis of a ratio 1/2 between the energy con- 
tent in CR and in thermal gas obtained with numerical 
simulations of the acceleration processes. All the more so, 
when account is taken of the fact that the value adopted 



for comparison is the one reported by Bowyer et al. (1999), 
which is a factor 2.7 lower than the one later presented in 
Bowyer et al. (2004). 



7. Discussion 

It is convenient to start by focusing the discussion on the 
central value Bq and the radial behavior of the magnetic 
field. The values of Bq in Table [3] span an interval from 
1.2 (Coma) to 4.6 fiG (A2319). In all four cases, the cor- 
responding pressure is well below the CR pressure. If the 
"true" magnetic field strength were substantially smaller 
or larger, the consequence for the SRE in the central re- 
gion would be, in the first case, that their contribution 
must be negligible (otherwise £o should be unacceptably 
large), in the second case that their contribution is even 
more likely to be important, since £o could take a value 
lower than 0.3, which is the reference value adopted above 
to approximately represente the "maximum" RP content. 
A comparison with the current estimates of Bq from FR 
measurements is therefore in order. 

None of the values of Bq in Table [3] seems to be in- 
consistent with typical FR estimates in the inner regions 
of clusters, but one must take into account (see Paper I 
and references therein for a brief presentation of this is- 
sue) the rather large uncertainties affecting the FR results, 
which are mainly related to the approximations adopted 
for the spectrum of the field fluctuations as a function of 
their scale length. In particular, the value for the Coma 
cluster is in fair agreement with the estimate obtained by 
Kim et al. (1990) using background sources, but not with 
the one obtained by Feretti et al. (1995), using one ex- 
tended source within the cluster, which could be as large 
as ~ 8 /iG, if a very small scale length for field rever- 
sals is adopted, of the order of 1 kpc. The value recently 
derived for A2255 (Govoni et al. 2006), using three ex- 
tended sources within the cluster and a more sophisti- 
cated approach than in Feretti et al. (1995), is 2.5 /j,G 
(the confidence level is not explicitely given in the pa- 
per), which is not far from the 4.0 /iG, and practically 
equal to Bq, in Table [31 For the two other clusters the "in- 
dividual" FR measurements available are not reliable. It 
must, in addition, be noticed that, in the frequency range 
around 1.4 GHz, the SRE spectral energy slope is about 
3.7 (for the choice made of s = 2.7), because of the radia- 
tive losses, and the radio emissivity scales as ij( 3 - 7 + 1 V 2 
= B 2 - 35 : thus, it is very sensitive to the B value, as also 
illustrated in Figs. [1] to 2J Therefore, so long as a very 
robust FR estimate for each individual cluster is lacking, 
a generic agreement is by no means a strong enough ar- 
gument for concluding that the SRE might always be an 
important, if not the dominant, contribution. Moreover, a 
very relevant uncertainty on the values of B for the three 
clusters other than Coma stems from the large uncertainty 
on the radio spectral slope, hence on their specific value of 
s and consequently on the RP pressure for a fixed value 
of £o (Eqs. [31 [H and [5]). At present, it can only be stated 
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that the above-mentioned conclusion cannot be excluded 
for the inner region of a cluster and its radio halo. 

The situation is plainly different at distances from the 
center that are much greater than r c th- I n the three ob- 
jects where is several times larger than the core radius, 
Figs. [U [21 and|3]show that the magnetic pressure, although 
it remains well below the CR pressure, tends to increase, 
a behavior that is physically embarrassing. Rather than 
sharply concluding that the SRE should therefore be re- 
garded as a minor component of the bulk RE, it is perhaps 
more sensible to consider the possibility that the contri- 
bution of the PRE might become dominant at large radii, 
thus providing a way to get around this problem. To qual- 
ify this possibility, one can, for instance, assume a con- 
tinuous injection of fresh PRE throughout the cluster, 
described by a source function Qep(lei' r ')i with a radial 
dependence of the form (8) adopted above for the accu- 
mulated RP and an energy dependence with the same in- 
dex s. The solution (14) would then give the equilibrium 
distribution, whose energy dependence is practically the 
same as that of the SRE. One of the results for the case 
of Coma is illustrated in Fig. \7\ Here the PRE and SRE 
contributions at the center are assumed to have equal 
weight (that is, the quantity £o was reduced to 0.15 to 
ensure that the SRE central density amounts to half the 
value derived previously, while the value of Bq obviously 
remains the same, 1.2 /iG). The rise in magnetic pres- 
sure present in Fig. [1] has gone away almost completely, 
except for the mild bump around 400 kpc, which is evi- 
dently the artificial consequence of r c n being about twice 
r c ,th (see Tables Q] and [2|) , together with the prescription 
t c ,rp — fc,th adopted in Sect. 3. In conclusion, to obtain 
an acceptable radial behavior of the magnetic field, one 
must resort to PRE at large distances from the center, 
but they need not necessarily be the dominant fraction 
within r c tht where the production of SRE is most effec- 
tive. 

Observationally, the case of the UV excess analyzed 
in Sect. 6 is very instructive. If the surface brightness of 
the IC emission by the electrons responsible for the ra- 
dio emission could be measured in the hard X-ray band, 
the same procedure could be followed to verify whether 
or not the amount of parent RP is compatible with the 
condition £ less than unity everywhere throughout the 
cluster. A positive answer would not, however, be con- 
clusive as to the secondary origin of the RE. This issue 
can only be addressed directly through measurements of 
the collisionally-produced GR. This brings one back to the 
key issue of the effective CR content of clusters. 

In this respect, a final important remark is due on the 
apparent dichotomy between those clusters endowed with 
a radio halo and those that are not. In the current liter- 
ature (Feretti, Burigana & Enfilin 2004), a correlation is 
emphasized between the existence of a radio halo and the 
evidence of ongoing merging phenomena. This correlation 
speaks in favor of a dominant population of PRE every- 
where in these clusters, which are being accelerated since 
at most one billion years. However, it must be stressed 
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Fig. 7. Upper panel: hypothetical mix of SRE, in the case 
£o = 0.15, with PRE of equal density at the center of the 
Coma Cluster. The density refers to electrons with -y e — 10 4 , 
the energy corresponding approximately to the synchrotron 
frequency v c =\A GHz when Z?o=1.2 \iG. Lower panel: the ra- 
dial behavior of the magnetic field required to reproduce the 
radio surface brightness. See text for details. 



that an attempt to estimate the strength of the magnetic 
field in those clusters devoid of a diffuse radio emission 
by combining FR measurements in an ensemble of them 
(Clarke, Kronberg & Bohringer 2001) resulted in a value 
of B in the interval 1 to 10 /iG, the same interval of the 
clusters with a radio halo. The results obtained in Sect. 
5 would then suggest that the amount of RP, which have 
accumulated in both types of clusters all along their ten 
billion years history, should be significantly lower than the 
amount corresponding to the reference value £o = 0.3. It is 
the opinion of the authors of the present paper that more 
sensitive measurements of a diffuse radio emission in clus- 
ters, where it has gone undetected so far (along with more 
accurate "individual" estimates of B using the FR) , might 
provide very stringent constraints on their CR content. In 
consideration of the hard fact that the detection of GR 
from more than just one cluster will be a difficult task even 
for the mission GLAST, further efforts should be made in 
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the radio window. In order to stimulate such efforts, some 
clusters, where a diffuse radio emission has gone unde- 
tected so far, and which look like promising targets, have 
been selected (Table [5J lower part). For reasons of ho- 
mogeneity with the ones discussed in this paper, none of 
these clusters shows the soft excess usually named "cool- 
ing flow" . The diffuse radio emission from these clusters, 
which should be due exclusively to SRE, should be typi- 
cally more concentrated than in the "radio halo" clusters, 
because of the dependence of the production rate from the 
product n t h X ripp. Thus, the 1.4 GHz flux (Table O col. 
3) was calculated (according to the prescriptions in Sect. 
3) for each of them and on the basis of its X-ray prop- 
erties, from within r = r c t \ l (Table col. 2), assuming 
the reference value 0.3 for £ , and a value equal to 5 /xG 
for J5 . Note that the predicted flux values scale linearly 
with £o and, in first approximation, with the square of 
Bq. In the same table, for comparison, the flux is reported 
from the corresponding region of the four clusters studied 
above. Despite the large uncertainties associated with the 
actual value of Bo, it is encouraging that the predictions 
are comparable, in order of magnitude, with the measured 
values of the four clusters. Furthermore, as reported in the 
same table (col. 4), interesting upper limits already exist 
for three of them. One case seems particularly significant, 
namely A119, where the prediction and the upper limits 
are comparable, while some FR measurements are avail- 
able (Clarke, Kronberg & Bohringer 2001; Feretti et al. 
1999) that indicate a Bq value of about 5 fj,G in this clus- 
ter (Murgia et al. 2004). Thus an improvement by a factor 
of a few in the radio upper limit could already place a sig- 
nificant constraint on £o- Also the case of A399, where the 
upper limit is about ten times lower than the prediction, 
could be very significant, but the authors are not aware 
of reliable FR measurements for this object. The last col- 
umn of Table [5] reports the predicted gamma-ray fluxes 
in the 0.1 — 10 GeV interval from within r = r c ^h'- from 
a comparison with the similarly defined fluxes of the four 
clusters in the upper part of the table, they appear to be 
hardly detectable, with the possible exception of A3667. 



8. Conclusions 

The issue of the CR content in galaxy clusters, which are 
endowed with a diffuse radio emission called radio halo, 
was discussed here in terms of the secondary products of 
collisions with the thermal gas, namely SRE and GR. 
This production is maximized when the density profile of 
the RP is identical in form to that of the thermal gas. The 
production rate can be normalized by adopting a value of 
the ratio between RP and thermal pressure, £oi constant 
throughout the cluster out to the radius Rh of the radio 
halo. This quantity represents the efficiency and output of 
current models of CR acceleration associated to merging 
processes in clusters. The value adopted here as a reference 
is £o = 0-3 to represent the case of maximum efficiency 
and cumulation over the cluster lifetime. A representative 



value s — 2.7 has been chosen for the energy spectral index 
of the RP. 

The equilibrium spectrum of the SRE is then con- 
strained to reproduce the observed radio emission, and 
both the central value Bq and the radial behavior of the 
magnetic field is then uniquely derived. In the four clusters 
treated, the values of Bq fall between 1.2 and 4.6 fj,G (re- 
duced to 0.8 to 3.2 fiG if scalar fluctuations on the order of 
< (SB) 2 > /B 2 — 1 are admitted): these values are con- 
sistent with estimates of the same quantity from FR mea- 
surements within clusters in general, but more detailed 
FR measurements are needed for each cluster individually 
before a sound conclusion might be drawn. Particularly in 
view of the high sensitivity of the synchrotron emissivity 
to the value of B. 

The radial behavior of B, instead, seems unnatural, 
in that it tends to increase with the radial distance from 
the center. Rather than definitely exclude a relevant con- 
tribution of the SRE to the radio emission also in the 
central parts of a cluster, this result suggests that PRE, 
accelerated "recently" , i.e. within their maximum radia- 
tive lifetime, must dominate at large radii. Figure [7] shows, 
in the case of the Coma cluster, how the mix of the two 
contributions alleviates the problem of the radial behav- 
ior of B, even if the contribution of the SRE within the 
thermal core radius is assumed to be as large as that of 
the PRE. 

The GR production rate in the same conditions de- 
scribed above immediately yields the expected luminos- 
ity and flux in the 0.1 to 10 GeV. A comparison with 
the sensitivity curves of the future missions AGILE and 
GLAST shows that, out of the four clusters considered 
here, only the Coma cluster could be very significantly 
detected by GLAST, even if £ were as small as 0.1. With 
£o = 0.3, Coma could be marginally detected by AGILE, 
while A2319 and A2163 could be marginally detected by 
GLAST. 

The proposal by Bowyer et al. (2004) that the 
spatially-resolved UV excess emission in the Coma 
Cluster, supposedly due to IC scattering by RE on the cos- 
mic microwave background, might originate from SRE, is 
discussed as an instructive example of what could be done 
if spatially-resolved measurements of the IC emission from 
the radioemitting electrons, expected in the hard X-rays, 
were available. The SRE hypothesis is rejected on the 
grounds of an estimate of £ that, at the center of the clus- 
ter, should be at least one order of magnitude greater than 
unity. 

Finally, the dichotomy between clusters with and with- 
out radio halos is presented as an opportunity to test the 
CR content in clusters in general. This follows from the 
fact that clusters without a diffuse radio emission, accord- 
ing to FR measurements (Clarke, Kronberg & Bohringer 
2001), share a magnetic field of similar strength with those 
endowed with a radio halo. The results show that, if £o 
were close to the reference value adopted in this paper, 
the SRE would give rise to a measurable radio emission. 
Therefore, any effort to improve the limits on diffuse ra- 
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Table 5. The four clusters studied in this paper (upper part) and a sample of six clusters whitout a radio halo (lower part) 



Cluster 


1"c,th 


Fi. 4 (r < r^th) 


i*i. 4 (r < r c ,th) 


F g {0.1 - 10 GeV,r < r c>tk ) 




(arcmin) 


(mJy) 


(upper limit, mJy) 


(cm -2 s _1 ) 


Coma 


10.5 


160 




6.8 x 10" y 


A2319 


1.6 


18 




2.7 x 10" 10 


A2255 


4.8 


25 




2.7 x 10~ 10 


A2163 


1.2 


18 




2.9 x 10~ 10 


A119 


5.3 


30 


< 47 


1.8 x 10 _iu 


A2063 


1.3 


6 




3.8 x 10" 11 


A1775 


1.8 


6 




4.1 x 10" 11 


A1413 


0.7 


8 




5.8 x 10" 11 


A399 


1.9 


22 


< 1.6 


1.5 x 10~ 10 


A3667 


7.9 


54 


< 209 


3.3 x 10~ 10 



The X-ray parameters adopted for the predictions in the lower part come from: Cirimele, Nesci & Trevese (1997) for A119, 
A2063, A1775 and A1413; Sakelliou & Ponman (2004) for A399; Vikhlinin, Markevitch & Murray (2001) for A3667. The flux 
upper limits in A119, A399 derive from Giovannini et al. (2006), the one in A3667 from Feretti (2004). 



dio emission from clusters in general, and the estimate of 
their magnetic field through FR measurements, will help 
to constrain their CR content before the hard task of the 
future gamma ray experiments can be accomplished. 
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